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ABSTRACT 

A Ground Winds Test was conducted by ^ the Marshall Space Flight 
Center, Unsteady Aerodynamics Laboratory on the IB/LC-39 vehicle in 
the fall of 1971, BAI provided on-site engineering at the Langley 
Research Center, Dynamics Wind Tunnel and later final reduction of * 
the tape recorded data. The Ground Winds Data Reduction System (GWDRS) , 
specifically designed for rapid and inexpensive data analysis was 
utilized to analyze these tapes for full scale, first and second mode 
bending moment or acceleration plots versus dynamic pressure or wind 
velocity. 

Select portions of the Skylab 2 tape data was analyzed statistically 
in the form of power spectral densities, autocorrelations, and cross - 
correlations to introduce a concept of using system response decay as 
a measure of linear system damping. The Engineering Notes written 
also discussed applying these techniques to nonlinear systems. Two 
small test instruments utilizing the latest LSI circuitry, namely 
Log-Decrement Dampometer and Auto correlator Dampometer, were designed 
and constructed for evaluation. 
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1.0 INTRODUCTION 

The on-site engineering suppoft provided at the Langley Research 
Center Wind Tunnel during the Satiim IB/LC-39 Ground Wind Tests and 
subsequent final reduction of the first mode strain gage data using the 
Ground Winds Data Reduction System were summarized in BAI Annual Report 
No. 605. This report also documented the Log-Decrement Dampometer 
designed and utilized during these tests. 

Final analysis for second mode full scale loads using the Ground 
Winds Data Reduction System on tape recorded second level strain gage 
signals was performed under contract Amendment (DCN 1-2-75-20055). 

Plots of dynamic lift, drag and static drag bending moments versus wind 
velocity or dynamic pressure were produced for 112 Data Points encompass- 
ing five fueled vehicle configurations. These final results were trans- 
mitted to the OOR in an original and one copy of BAX Data Report 607-A. 

Final analysis for first mode lift, drag, and resultant acceleration 
versus dynamic pressure utilizing the GWDRS on the tape recorded x and y 
accelerometer data at Body Station 121 was continued under Contract 
Amendment (DCN 1-3-75-30019). An original and two copies of BAI Data 
Report 608, Volumes I through VI were mailed to the COR as the data 
reduction was completed. 


(l.D 


v 



BAGANOFF ASSOCIATES, INC. ' 

Report No. : 610 

Date: January 18, 1974 


Statistical analysis in the form of narrowband power spectral 
densities was conducted on a total of 80 Data Points - tape track 
combinations to further define the first modal frequencies. These 
results in the form of digital tabulations and plots were mailed to ' 
the COR for use in follow- 70 wind tunnel tests conducted in the labora- 
tory. 

During this contract, RAI perfoimed introductory studies in the use 
of computer technology for statistical analysis of multi-modal aero- 
elastic systems. See Engineering Notes listed in Section (3.2). This 
effort culminated in the design and construction of an Auto correlator 
Dampometer for calculating the system damping ratio under random excita- 
tion conditions. A computer simulation program was first written to 
define the instrument’s design parameters before construction began by 
applying the range of damping ratios encountered in typical wind tunnel 
tests. 

With the advent of micro-processors, a continuous design effort is 
necessary between the aerodynamicist and the computer technologist in 
order that the former may avail himself of new measuring techniques 
not generally available. Results on the evaluation of the Autocorrelator 
Danpometer on Saturn IB/LC-39 Ground Winds Test data are contained in 
Section (4.0). 


( 1 . 2 ) 
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2,0 DATA REDUCTION 

2.1) Skylab 2, Second Mode Results - The second level strain gage signals 
were reduced using the Ground Winds Data Reduction System (GWDRS) for full 
scale dynamic lift and drag bending moments versus wind velocity and 
static drag bending moments versus dynamic pressure. The following 112 
Data Points were plotted and presented in BA.I Data Report 60 7 -A: 


74 through 82, 95 through 110: Empty Weight (Secondary Scaling), 

161 through 167, 181 through 191: Intermediate Weight (Secondary 



Scaling) RP1 in SIB State, Space- 


craft Fueled. 

270 through 322: 

Completely Fueld (Secondary Scaling). 

328 through 349: 

Completely Fueled (Primary Scaling). 

390, 402: 

Empty Weight (Primary Scaling) . 


A representative plot of full scale Dynamic Lift Bending Moment versus 
Wind Velocity is shown in Figure (2.1). The GWDRS also produced the full 
scale Dynamic Drag Bending Moment versus Wind Velocity plot in Figure (2,2) 
and the full scale Static Drag Bending Moment versus Dynamic Pressure plot 
in Figure (2.3). 




SKYLAB 2 GROUND WINDS TEST DATA 



• WIND VELOCITY AT STATION 1500 (KNOTS) 

CONFIGURATION: (XUPLETELY FUELED (SEOONDAKY SCALING) " £- . 

SEOCND MODE DYNAMIC LIFT BENDING MOMENT AT STATION 1077.27 0 

WIND TUNNEL A2IMUIH ANGLE 105° 

DAMPING 1.0 I, DATA POINT 285 
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WIND VELOCITY AT STATION 1500 [KNOTS) 
CONFIGURATION: COMPLETELY FUELED (SECONDARY SCALING) 

SECOND MODE DYNAMIC DRAG BENDING MOMENT AT STATION 1077.27 • 

WIND TUNNEL A2BJUTH ANGLE 105° 

DAMPING l.o I, DATA POINT 285 
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^ DYNAMIC PRESSURE (PSF) 

(6 

^ CONFIGURATION: COMPLETELY F UELED (SECONDARY SCALING) 

STATIC DRAG BENDING MOMENT AT STATICS! 1077.27 
W WIND TUNNEL AZIMUTH ANGLE 105° • 

DAMPING 1.0 I, DATA POINT 285 
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2 t 2) Skylab 2, First Mode Accelerometer Data - The x and y accelerometer 
signals derived at Bod/ Station 121 were reduced on the GWDRS to Resul- 
tant, Lift and Drag Acceleration versus Dynamic Pressure for Data Points 
70 through 157, Empty Weight (Secondary Scaling) Configuration. Repre- 
sentative plots are presented in Figures (2.4), (2.5), and (2.6). These 
results are contained in BAI Data Report 608, Volumes 1 through VI. 


(2.5) 
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DYNAMIC PRESSURE (PSF) 

CONFIGURATION: EMPTY WEIGHT (SECONDARY SCALING) 

RESULTANT ACCELERATION AT STATION 121 
WIND TUNNEL AZMJIH ANGLE o° 

HOPING 5.0 %, DATA POINT 70 
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WIND TUNNEL AZIMUTH ANGLE 0° 
DAMPING 5.0 %, DATA POINT 70 
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3.0 STATISTICAL ANALYSIS 

3.1) Skylab 2, Power Spectral Densities - The Narrowband Analog, Power 
Spectral Density System was utilized to produce results with a small 
statistical error. Since a bandwidth of 0.25 Hz was required in order 
to define the first mode frequencies, an averaging time of 60 seconds 
was utilized. In all eighty analyses were conducted with a PSD analysis 
performed at four intervals A, B, C, and D during a sweep Data Point. A 
typical digital tabulation and corresponding graph are presented in 
Figures (3.1) and (3.2). 

3.2) Skylab 2, Introductory System Damping Studies - Three Engineering 
Notes were written during this period presenting the results of statis- 
tical analysis in terms of autocorrelations and cross- correlations of 
Skylab 2 data for system damping information. It is firmly believed 
that statistical techniqes, such as these, will help to formulate an 
engineering model for predicting an understanding of self- excited system 
loads . 

1) ’’The Behavior of the Cross -Correlation Function For Viscous 
Type Damping 11 , January 14, 1972. 

2) " Correlation Analysis of Self- Excited Systems”, Feb. 28, 1973, 

3) "Descriptive Model, Ground Wind Loads,” March 27, 1973. 

(3.1) 
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BAGANOFF ASSOCIATES, INC 
SKYLAD 2 GROUND WINDS 
PSD FIRST MODE. 

DATA PT.104 
TAPE TR. 2 

, RUN SEC. C 

AVE- TIME 60.4 Sec. 
B «= .25 HZ 

FILTER CTR. ... PSD 

FREQ. CHZ ) CV)SQ/HZ 


1 

. 1 500E+0a 

. 1066E+03 


1 0 6 6 F + 0 3 

. 1-581E-05 

. 158 IE- 05 

e 

♦ 1535E+02 

. 1075 E +0 3 


1 101F.+0 3 

. 1944E-05 

. 1 944F.-0 5 

3 

. 1 5 50 E+03 

• 1 080E+03 


1 1 8 2 E 4 0 3 

.218712-05 

.R187E-05 

A 

• 1575E+03 . 

• . 1076 E + 0 3 


1 1 3 5F. + 03 

. 1974E-05 

. 1974E-05 

5 

. 1 600E+02 

. 1 1 05E403 


1 152E+03 

•3866E-05 

* 3 8 68 E- 05 

6 

. 1 625E+02 

. 1105E+03 


1 165E+03 

•3 8 61 E -05 

•3861 E -05 

7 

. 1650E+03 

.1 1 i 9E + 03 


1 1 78E+03 

• 5408E-0 5 . 

• 5408E-0 5 

8 

. 1 675E+03 ' 

. 1 1 1 8E + 03 


1 187E+03 

• 4 56 7E-05 

• 4 5 6 7E - 0 5 

9 

. 1700E’+02 

. 1 127E+03 


1 1 96E+03 

.643815-0 5 

• 6438E-05 

10 

. 172SE402 

..1 1 33E + 03 


1206E+03 

. 7424E-05 

. 7 4 2 4 E - 0 5 

l i 

. 1750E+02 

.1 1 72E+03 


1 222E+03 

. 1 793E-04 

. 1 79 3 E- 04 

12 

. 1775E+02 

. * 1 1 72E + 03 


1234E+03 

. 1831E-04 

• 1 8 3 1 E - 0 4 

13 

. 1800E+02 ■ 

. • 1 1 86E+03 


1 246E+03 

• 25 1 5E-04 

. 2515E-04 

14 

. 182SE+02 

. 1046E+O3 


1 276E+03 

• 9 9 9 1 E - 0 4 

.999 IE- 04 

15 

. 1B50E+02 

. 1S81E+03 


1309E+03 

.221 IE-03 

.281 IE-03 

16 

. 1875E+02 

■ 1 296E+03 


1 333E+03 

•3175 E -03 

.31 75E-03 

17 

. 1 900E+02 

. 1267E+03 


1342E+03 

. 1599E-03 

. 1 599E-03 

18 

• 1935E+03 

•1831 E+03 


1 3 A 5 E + 0 3 

• 7030 E -04 

. 70 30 E- 04 

19 

. 1 950E+02 

. 1 196E+03 


1346E+03 

. .3165E-04 

• 3 1 65E-04 

20 

. 1975E+02 

.1 183E403 


1 34 7E! 1 03 

•2346 E -04 

• 2 3 4 6 H> 0 4 

21 

• 2000E+02 

.1 181 E+03 


1348E+03 

• 2232E-04 

• 2232E-04 


Figure £3*1) 
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The promising test instrument described in Section (4,0) incor- 
porating ROM’s, mini -processors, and RAM’s represents a look into the 
future of wind tunnel testing. 
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4 , 0 ALTOCORRELATOR PAMPOMETER 

The system dancing associated with a wind tunnel model of a launch 
vehicle subjected to increasing ground winds is an important parameter 
from the standpoint of monitoring the integrity of the model, and learn- 
ing about the self- excited flow process. BAI has been in the forefront 
of developing methodology for using statistical analysis in the form of 
autocorrelation functions, etc. for the gaining of insight into the. 
closed- loop system composed of the aerodynamic forces and the launch ve- 
hicle response. See Engineering Notes listed in Section (3.2). As is 
the practice in servomechanisms, measurements can be made on the closed- 
loop system and characteristics of these measurable quantities attributed 
to the physical parameters of the various elements in the closed- loop, in 
this case the aerodynamic force and the launch vehicle response. 

. The Auto correlator Dampomoter described herein was developed under 
this Contract and represents a small and versatile test instrument that 
can measure the system damping on-line in the wind tunnel with the wind 
blowing. This instrument is designed such that when the occasional high 
peaks in the random response exceed the voltage threshold at the input 
to the instrument, an input gate enables peak detection circuits to 


(4.1) 
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apply to memory the peak values for the succeeding 32 vibration cycles. 

These 32 peak values would, in theory, conform to a cosine exponential 
decay except that the subsequent random inputs to the system cause these 
peaks to not conform to the free vibration envelope. The electronics in 
the heart of the instrument multiplies together the peak values and sums 
these products for succeedingly large delay times tau in the sense of the 
ordinary autocorrelation function, and the output becomes the systematic 
impulse response for the system with the effects of the random fluctuation 
of the input eliminated. This instrument has the advantages that: (1) it 

operates with peak values which maximize the signal- to -noise ratio and are 
more easily detected electronically. And (2) the large energy in the 
peaks is primarily that of the desired vibration mode .in a multi-modal 
system. 

An analytical model for a ground winds vehicle falling into the cate- 
gory of self-excited systems is practically non-existent, thus each time a 
now launch vehicle is designed, a wind tunnel test must be conducted in 
order to measure experimental loads for various launch configurations and 

wind velocities and orientations. A predictive model does not exist because 
empirical data in the form of statistical averages are a rarity since 
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reliable test instruments for making on-line measurements on a random 
process are just coming into being. 

Off-line, large digital computers are capable of performing the 
matrix manipulations required in any given statistical analysis once the 
analog signals have been preconditioned and digitized. Typically , the 
large computer is required to handle large quantities of data points , and 
the results are available long after the wind tunnel tests have been 
couple ted. The memory capacity of even a disc drive on a mini -computer 
is soon exceeded if the incoming analog signal is sampled continuously 
according to the Nyquist criteria. The Autocorrelator Dampometer may be 
thought of as consisting of two units, first a data compressor to pass only 
the segments of the data required, and second a matrix manipulator to per- 
form the desired statistical analysis*. 

The standard procedure in the wind tunnel for obtaining structural 
damping data is to excite the model sinusoidally with on electromechanical 
exciter, and then record the free vibration decay on an oscillograph. 

The desired information is obtained from the log- decrement of the decay 
envelope. The problems with this method are: (1) expensive wind tun- 

nel time is lost while hand calculations are made,, and (2) in multi-modal 


(4.3) 
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systems the decay may pertain to the desired mode for the first few 
vibrations, and the later smaller vibration cycles pertain to another 
mode. 


Under Contract NAS8-26703, DCN 1-1-75-10052 (IF), an instrument call- 
ed a Log-Decrement Dampometer was designed and constructed that electroni- 
cally sampled the peak values of the free vibration signal and computed 
the log of the ratio of successive peak values, and displayed in decimal 
notation the average structural damping ratio. This instrument, which 
occupies approximately 1,3 cubic feet, was utilized in the Ground Wind 
Tests conducted at LRC on the Saturn IB/LC-39 vehicle in the fall of 1971. 
This Log-Decrement Dampometer reported in BAI Report No. 605 was found to 
have its limitation in the wind tunnel environment because, (1) the unit 
had to be timed with the onset of the free vibration, (2) a multi-modal 
decay produced an average damping ratio which differed slightly from the 
desired average, and the supposedly free vibration always contained a 
random noise background due to compressor motor inputs. The lessons 
learned using this instrument helped develop added impetus to the al- 
ready held concept for the Auto correlator Dampometer. In Figure (4.1) 
the Log-Decrement Dampometer is the unit shown on top in the picture, 
and the more recently constructed Autocorrelator Dampometer is 
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Figure (4.1) 


(4.5) 
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AUTOCORRELATOR DAMPOMETER EVALUATION Date: January 18. 1974 

Dwell Data Point 209 Q = 50 PSF 
$ = 120 ° Damping Ratio = .421 


Input: Time History of 

Base x Strain Gage 




Output: Autocorrelation 

r _ . Function Repeated 
< I /- 0 
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Dwell Data Point 214 Q = 120 PSF 
<t> = 120° Damping Ratio 


Input: Time History of 



(4.21) 


- 1.17% 

Output : Autocorrelation 






i Figure (4.4) 
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the unit shown on the bottom. 

Under the current Contract Amendment, the Autocorrelator Dampometer 
was designed and constructed which is capable of operating in the wind 
tunnel while the wind is blowing and thus produce the total system damp- 
ing ratio. The next section discusses the definition of the auto- 
correlation function in series fom and gives results for a computer 
simulation encompassing damping ratios from 0.51 to 5% as encountered in 
the ground wind tests. The simulation showed that the Autocorrelator 
Dampometer was capable of producing a more ideal linear system, free 
vibration response as an output even when the input contained two and three 
slope changes, and also noise modulation. This study determined the 
optimum number of digital registers required in the unit to be constructed. 


( 


(4.6) 
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4.1) Simulation of the Autocorrelation Function - A computer program 
was first written to evaluate the effectiveness of the autocorrelation 
function in determining percent critical damping for the range of values 
experienced in the Ground Wind Tests. Accurate solutions were obtained 
for structural damping ratios from 0.5% to 5% and three slope changes 
during a case single simulated decay. 

4.1.1] Autocorrelation Function - The autocorrelation function for 
stationary signals may be written as: 


X 1 

R ffCr)“ 4-J fCt)fCt-T)dt (4.1) 

Where: R is the autocorrelation function 

T is the integrating (averaging) time 
f is the waveform of the time function 
t is the delay or time shift 

If f(t) is a sampled signal, then Equation (1), (provided N is relatively 
large) can be approximated by: 

R ffO< T s ) = f(nT s )£(nT s - kT s ) (4 . 2) 

Where: N is the total number of sample products 

k is the total number of sample period shifts 
T^ is the sample period 

For special application to vibration decay signals, the waveform' 
is an exponentially decaying sinewave. Only the input waveform peaks 
need to be sampled in order to retain the critical damping information. 

The output waveform for the autocorrelation function is likewise an exponentially 

(4.7) V 
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Where: represents the pea k reading o£ the nth sample. 

An exponentially damped sinewave having a decay rate corresponding 
to 2% critical damping was selected as the most likely rate to be encountered 
in wind tunnel tests. Peak readings for 26 cycles of this waveform, serving 
as input to the computer program, are tabulated in Table (4,2) and plotted 
in Figure (4.2a). Decay rates corresponding to 0.5% and 5.0% damping were 
selected as limiting cases to be found in actual tests and peak values 
for these input waveforms are also found in Table (4.1). To. simulate a 
more realistic waveform as experienced in actual tests, an additional 
case was generated where the effective damping was allowed to change 

after eight cycles from 2% to 1% damping as shown on the graph on Figure (4.2b). 
Still another case was generated by allowing the decay to change again 

after 12 cycles back to 2% damping. (4.8) 




Ptm IK. 
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Peak values for this waveform are tabulated in Table (4.1) and depicted 
in Figure (4.2c). A third test case was made by applying a 2 % modulation 
to the first decay rate. These values are also tabulated in Table (4.1). 


Table (4.1) iNFUi FEAKT VALUfccrTUKTEol 


Sample No. 

An 

•0*5% Damping 

2 % Damping 

5 % Damping 

2 % to 1 % 
Damping 

■ 

2 % Modulation 

0 

1 — 

1,000 

1.00 

1.000 

1.00 

1.00 

1.00 

1 

.965 

.88 

.720 

.88 

.88 

,8976 

2 

.940 

, .77 

.540 

.77 

.77 

.7546 

3 

.915 

.68 

• .394 

.68 

.68 

.6936 

4 

.880 

.60 

.286 

.60 

.60 

.588 

5 

.855 

.53 

.210 

.53 ‘ 

.53 

.541 

6 

.825 

.47 

.152 

.47 

.47 

.46 

7 

.800 

.41 

.110 

.41 . 

.41 

.418 

8 

. .770 

.365 

.079 

.365 ' 

.365 

.37 

9 

.745 

.32 

.058 

.34 

.34 

.32 

10 

.720 

.285 

.043 

.323 

.323 

,28 

11 

.700 

.25 

.031 

.305 * 

.305 

.25 

12 

. .685 

.222 


.287 

.287 

.21 

13 

.665 

.198 


.368 

.210 

,198 

14 

1 .645 

.174 


.25 

.197 

.174 

' 15 

.625 

.150 


.235 

.184 

.150 

16 

' .605 

.134 


.22 

JL71 

.134 

17 

.587 

.118 


.205 

.160 

.118 

18 

.568 

.104 


.192 

.151 

.104 

19 

.550 

.092 


.180 

.142 

.092 

20 

.532 

.081 





21 

.515 

.071 





22 

.500 

.063 





23 

.487 

' .055 





24 

.470 




, 


25 

.458 






26 

.443 






27 

.428 






28 

.415 






29 

.402 






30 

.390 






31 

.377 

A ‘ 





* 32 

.364 





* 


| Computational Errors For Nominal 0.5%, 2.0% and_5% Dajfpirig - 

When the various autocorrelation points are computed, the last few 

(4.10) 
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products have zero values because of the time shift of the A^ samples. 

For example, in computing R(T]_) the last multiplication in the summation 
would be Ai 9 <A 2 q. Since A 20 has a zero argument this product will equal zero. 
An algorithm was derived experimentally that suggested the first sample 
be used over again for A 201 appropriately weighted by a constant, i.e. 

^20 “ Therefore, a new set of equations would result for R(t) and 

as shown below: 

R Ct 0 ) =20 n l 0 V** (4.8) 

R C^i) n*0 t A n 0 A n +l + Ai9«koA 0 ) (4.9) 

R(T 2 ) “2^ Jo (V Vt-2 +A 18' k o A o +A 19 k l A l) (4.10) 

RCt 3 ) Iff Jo (VV 3 +A 17 ' k oV A 18 ,k lV A U k W (4.11) 

R ( t 4 ^ n £ 0 ( A n' A n+4 +A 16‘ k o A o +A 17' lc l A l +A 18' k 2 A 2 +A 19* k 3 A 3) (4. 


The values of kQ through k 3 are calculated from the following equations: 


K = 


A20 


K 


0 ~^o 

A 21 

X “T 


K. 


^22 

A 2 


A« - 
K = -23. 

3 A 3 


(4.13) 

(4.14) 

(4.15) 

( 4 . 16 ) 


Using the input peak values for. Test Case 2, R(t q ) through R(t, ) were 

(4.11) 
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computed using both sets of Equations (4.3) through (4.7), and the 
Modified Equations (4.8) through (4.12). Values of R(t) f s are listed in 
Table (4.2) below: 

Table (4.2) 

R(t) Computations for 2.0$ Nominal Damping Data and 20 Memory 
Locations and Extended Memory 

R(t) R(t) 

20 Locations .. Extended; Data 



.2216 

.2216 

\ 

.1948 

, .1952 

■ 

.1712 

.1719 

* 

.1506 

.1518 

■ * 

.1324 

.1340 

' Computed Danping 

(2.00$) 

(2.00%) 


Both sets of computations: (1) summations containing 20, 19, 18 -t- 

products and (2) summations containing a constant 20 products, produced the same 

answer for the computed damping of 2.00$. The implementation using 20 fixed 

locations and no natural extension of the input data gave identically 

2.00$ damping because the quantities of A„_ — were small to begin 

17 20 

with. 

The extended data solution in Table (4.3) gives the almost perfect answer 
of 0.49$, meaning that 36 memory locations would be the most desirable. 

Thirty- two memory locations produce the next best result, of 0.64 $ 

damping and the remaining case 1.0$. Tentatively, an autocorrelator incor- 

(4.12) 
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porating 32 memory locations might be the best compromise if the results 
for 5%. critical damping are satisfactory. 

Table (4.3) 

R(t) Computations for 0.5% Nominal Damping Data . 
and 20 and 32 Memory Locations and Extended Memory 



R(t) 

20 Locations 

R(t) 

32 Locations 

R(t) 

32 Locations, 

T o 

.5845 

.4426 

.4426 

T 1 

.5517 

.4246 

.4289 

T 2 

' .5196 

.4070 

.4156 

T 3 

.4877 

.3898 

.4027 

T 4 

.4561 

.3728 

.3900 


Confuted Damping 
(1.0%) 

(.64%) 

(.49%) 


Since for 5.0$ clamping the A n r s decrease rapidly to small magnitudes 
20 memory locations produced a very satisfactory answer 4.98%, and 32 


data locations would not improve it. 



Table (4.4) 

R(t) Computations for 5.0% Nominal Damping Data "i. 
and 20 Memory Locations 

R(t) 

20 Locations 
.1069 
.0780 
.0576 

’ .0419 V 


(4.13) 
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4.1.4) Computational Errors for Changes in the Decay Slope - 

Peak magnitudes for this evaluation are tabulated under heading 
2 % to 1 % Damping, in Table (4.1). Using the unmodified Equations (4.3) 
through (4.7), (no extension of the input data) the resulting R(t)*s are 
listed in Table (4,5) 

Table (4.5) 

R(t) Computations for Nominal 2% to 1% 

Damping Slope Change 

R(t) 

20 Memory Locations 

t 0 .2367 

Tjl . 2096 

t 2 .1857 

t 3 .1651 

t 4 .1470 

Computed Damping 
' (1.971) 

The computation of structural damping for the primary mode of interest 

in a model, 1.97% is reasonably accurate with the contribution from the un- ■ 
wanted mode, 1% damping , practically excluded. 

4.1.5) Computational Errors For Decay Slope of Z% to 1% arid Back to 2% - 
The peak input values for this ‘evaluation are tabulated under the 

heading "3 Slopes”, in Table (4.1). Using Equations (4.3) through (4.7), 
the following solutions for R(x) were computed. 


(4.14) 
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Table (4.6) 

R(t) Computations for Decay Slope of 21 to 1% and Back to 2 % 

R(t) 

20 Memory Locations 
tq ..2301 

T 1 .2030 

t 2 .1793 

t 3 .1588 

,1409 

Computed Damping 
(1.97%) 

, '• i 

The computed critical damping of 1.97% is a good approximation tovthe 
nominal damping of 2.0% representing the damping of the vibration mode 
of interest in the structure. 

4.1,6) Computational Errors For Input Data Composed of 2 % Nominal 
Damping Plus 2 % Noise Modulation - 

For this simulation, the' values under the heading >'' 2 % Damping 1 ' in Table (4.1) 
were alternately increased and decreased by 2%. The resultant input data 
are entered in the last column of Table (4.1) and the answers for R(t) are 
listed in Table (4.7) following. 


(4.15) 
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Table (4.7) 

R(t) Computations for 2.0% Nominal Damping 
Plus 2% Noise Modulation 

20 Memory Locations 
t q .2227 

^ .1958 

t .1713 

2 • 
t 3 .,1515 

t * 1326 

4 Computed Damping 

( 2 . 02 %) 

As shown, the autocorrelation process reduced the 2% noise modulation of 
the input data to 1% modulation of its output. The S/N . improvement was 
probably better than this, but the slope on the semi -log paper could not 
be read more accurately. 

4.1.7) Autocorrelator Design Detected By Structural Damping Accuracies - 
From the computer simulated tests, it was demonstrated that the 
autocorrelation function is an effective means of reducing the noise in 
a vibration decay signal and at the same time preserving the structural 
damping information. It is also concluded that 32 data storage cells 
in the autocorrelator were sufficient to produce acceptable results for 
input signals corresponding to 0.5% to 5% critical damping. 


(4.16) 
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4. 2) Autocorrelator Daripometers Lbgic 'Diagrams - As shown in Figure 
(5.3), the unit contains a trigger level adjustment and Peak Detector 
at its input. The trigger level adjustment allows the operator to se- 
lect the magnitude of the incoming strain gage signal at which, the 
Digitizer Control passes on an enabling trigger signal to the Analog- 
to- Digital Converter. • The Peak Detector only presents an output when 
an incoming signal peak has occurred. The Analog -to -Digital Converter 
converts the following 32 peak analog voltages to digital words, and 
stores them in the Random Access Memory (RAM) . The Read Only Memory 
(ROM) represents a mini -programmer for manipulating the words in memory 
and calculating the autocorrelation function. Complete documentation 
of the Auto correlator Dampoineter is contained in Figure (5.2) in the form 
of Logic Diagrams, Circuit Diagrams, and Wiring Diagrams. 

4. 3j) Skylab 2 System Damping Results - Data Points 209 and 214 from the 
wind tunnel series represented dwell conditions of 50 and 120 PSF res- 
pectively. These two runs at an azimuth angle of 120° produced re- 
latively large vehicle responses and these were analyzed previously for 
damping information by both analog and digital computer techniques. See 
Engineering Notes. Data Point 209 related to a tunnel velocity giving 
the critical Strouhal frequency and a system damping ratio of 0.42% while 


(4.17) 


i 



RAGANOFF ASSOCIATES, INC. 


Report No. : 610 

Date; January 18, 1974 

Data Point 214 gave a system damping ratio of 1.17%. 

Figure (4,3) shows a time history of the strain gage signal 
applied to the Autocorrelator Dampometer and the autocorrelation function 
at its output for three different segments of Data Point .209, In each case, 
a pencil mark notes the extreme peak where the Analog- to -Digital Converter 
first began operation. The 32 subsequent peaks were then digitized. The 
right-hand time histories show the autocorrelator function repeated on a 
very time contracted scale. This repeating action is necessary for refresh- 
ing the image of the autocorrelation function on an oscilloscope for opera- 
tor viewing. 

The first time history segment (1-D) presented in Figure (4,3) for 
Data Point 209 depicts 32 successively decreasing peak values beginning 
with the pencil mark. The Autocorrelator Dampometer output supports these 
observations of the input by portraying a rapidly decreasing stairstep func- 
tion corresponding to t s 0 and the first seven time lags. Time segment 
(1-E) depicts rapidly decreasing peak values for a time, and then the peak 
values remain at constant amplitude. The Autocorrelation output in (IE) 
can be seen to duplicate these findings . In time segment (1-F) , it can 
be seen that the peak values decrease for a time, and then begin to in- 
crease before 32 peaks have been accounted for. The autocorrelation out- 


(4.18) 
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Dwell Data Point 209 Q = 50 PSF 

<j> - 120° 1 Dancing Ratio » .42% 


Input: Time History of 

Base x Strain Gage* "" 



Output : Autocorrelation 
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in (IF) depicts a rapid decay more in agreement with the decay, in the 
first part of the time history. Linear response systems produce expo- 
nential decays which can be related to system damping ratios. The Auto- 
Correlator output also applies to nonlinear response systems, and this 
instrument may present a method of studying them. 

The first time history segment (1-B) in Figure (4,4) for Data Point 
214 shows a very rapid drop in the peak values followed by a large number 
of peak values of constant amplitude. The corresponding autocorrelation 
function exhibits approximately constant amplitude peaks. Segment (1-C) 
shows peaks that decay for a greater number of cycles but that once again 
flatten out. The cycles in the autocorrelation function (1-C) correspond- 
ing to the zero and first seven time lags tau supports the time history by 
showing a more rapid envelope decay. The time history for Segment (1-D) 
depicts a consistent decay for 32 cycles and the corresponding autocorrela- 
tion function (1-D) supports this observation with a rapid envelope decay. 

Figure (4.5) shows an expended time scale for the autocorrelator out- 
puts already shown in Figure (4.4) for Data Point 214. The departure of 
the seven stairstep functions away from the ideal exponential decay asso- 
ciated with a linear system would serve as an excellent instrument for 
studying nonlinearities in the flow process. This is the ' first time 
that such an instrument has been designed and constructed. 


(4.20) 
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Dwell Data Point 214 Q = 120 PSF 
= 120° Damping Ratio = 1.17% _ , 

Input: Time History of Output: Autocorrelation 
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The trigger threshold can be set so that in a normal five minutes 
wind tunnel run, whether a sweep or dwell condition, approximately 75 
autocorrelation functions would be calculated. A boxcar or exponential 
decay averager type could be constructed to ensemble average the variable 
autocorrelator functions for a time period consistent with the stationa- 
rity of the flow process. The construction of an autocorrelation ave- 
rager remains to be recommended as future work . The ensenbled auto- 
correlation function would then produce a hard number for the system damp- 
ing ratio in real-time for monitoring on an oscilloscope. 

The envelope decay or log -decrement of the ensembled autocorrelation 
function should produce an unbiased estimate of the true system damping 
ratio. This is based on the intuitive reasoning that the turbulent boun- 
dary layer which inparts a random force to the model at times constructively 
aids the vehicle motion, and at other times destructively interferes with 
this motion. There should be no preference for either constructive or 
destructive interference by the aerodynamic forces involved. 

The Auto correlator Dampometer output can be directly applied to 
the Log -Ratio Dampometer input as previously discussed for sinusoidal step 
function applications occurring in the wind tunnel and laboratory. In 
this case the latter instrument averages the log -ratio between successive 


(4.23) 
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steps of the autocorrelation output. The method of interconnecting 
those two instruments is shown in Figure (5.1). The Operating Instruc- 
tions are contained in Figure (5.2). The Computer Simulation Study of 
the Autocorrelator Dampometer presented in Section (4.2) indicates the 
power of utilitizing this instrument also for multi-modal analysis appli- 
cations under sinusoidal step inputs. 

Further progress in the de\ r elopment of these techniques requires . 
expertise on the part of the Unsteady Aerodynamics Laboratory at MSFC. 
and computer expertise as exists at BAI. 
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1. Connect cables as shown in figure 1. 

2. Remove top of Dampometer and do the following: 

A. Remove card 1 and replace it with card 12. 

B. Put switch on card 2 in A position. 

C. Put switch on card 9 in A position. 

D. Replace top of cabinet. 

3. Turn power on; both- units. 

SET-UP 

4. Set the front panel controls as follows: 

\ 

Correlator 

A. OPERATE- CAL switch in CAL. 

B. Depress RESET button. 

C. NEGATIVE PEAK off. 

D. Adjust input level to 1000. 

Dampometer 

A. LIMIT switch to 8. 

B . Dep re s s RESET ■ but ton . 

C. Adjust level adjust control for a reading 
of 1.500 on the Dampometer meter. 

5. Put Correlator OPERATE -CAL switch in OPERATE position. 

6. Set up for a solution is now ready. Set the- START and RESET 
buttons in the following order only: 

A. Correlator RESET 

B. Dampometer RESET 

C. ' Dampometer START 

D. Correlator RESET 

7. For high % damping where only a few cycles are available (above 4%) 
put the NEGATIVE PEAKS switch in the ON position. 

8. Adjust Correlator meter level to about 1000. 

9. Proceed as outlined from step 6. 

NOTES: 

1. A blank meter or a negative answer indicate ^.signal level 
too high to Dampometer. 

2. The CAL position of the OPERATE-CAL switch will put out 
a 5V signal on the Correlator output. 

The resolution when using negative peaks is 1/2 true solution. 
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Figure (5.16) 
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